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Abstract 


We present type IIB supergravity solutions which are expected to be dual to certain 
Lifshitz-like fixed points with anisotropic scale invariance. They are expected to describe 
a class of D3-D7 systems and their finite temperature generalizations are straightforward. 
We show that there exist solutions that interpolate between these anisotropic solutions in 
the IR and the standard AdS; solutions in the UV. This predicts anisotropic RG flows 
from familiar isotropic fixed points to anisotropic ones. In our case, these RG flows are 
triggered by a non-zero theta-angle in Yang-Mills theories that linearly depends on one of the 
spatial coordinates. We study the perturbations around these backgrounds and discuss the 
possibility of instability. We also holographically compute their thermal entropies, viscosities, 
and entanglement entropies. 
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1 Introduction 


The AdS/CFT correspondence has provided a very powerful and successful paradigm to analyze 
relativistic and isotropic fixed points in various quantum field theories [4]. On the field 
theory side, they are described by (d+ 1)-dimensional conformal field theories and are invariant 
under the homogeneous scaling transformation (t, £1, 2%2,...,2%q) > (At, Av, Av2,...,ATq). On 


the gravity side, they are equivalently described by gravity on a (d + 2)-dimensional AdS space 


d 
dr? 
2 2 2 X 2 


It is natural to try to extend the AdS/CFT correspondence to a holography for the following 


anisotropic spacetime 


p d 2 

dr 

2 2 2 2 2 2 
i=1 j=p+1 

where 0 < p < d — 1, and the parameter z(# 1) measures the degree of Lorentz symmetry 

violation and anisotropy. Since the metric (1.2) is invariant under the scaling (t, £i, yj, r) — 

(A7t, ATi, AYj, $), we expect that on the field theory side it is dual to a fixed point which is 


invariant under the scaling transformation 
(t, Ti, Yj) = (Mt, Ati, Ay;). (1.3) 


Notice also that by a coordinate redefinition r” = p, we can rewrite the metric (1.2) into another 


illuminating form (after rescaling (t, x;, y;) accordingly) 


P d 2 
d 
ds? = p* (a +5 aa?) $p X d+ 1 (1.4) 
i=1 j=p+1 p 


Thus we can equally argue that the dual background is invariant under an anisotropic scaling 
transformation (t, £i, Yj, P) > (At, Axi, Nü £), where the y; directions are responsible for the 
Lorentz symmetry violation and anisotropy. 

In general, fixed points with the anisotropic scaling property are called Lifshitz(-like) 
fixed ponid (see e.g. the textbook [7] for a brief review). This generalization of AdS/CFT 
correspondence to Lifshitz-like fixed points was first proposed and analyzed by Kachru, 
Liu and Mulligan in the particular case of p = 0. The simplest case with p = 0 represents 


non-relativistic fixed points with dynamical critical exponent z, which appear in many examples 


* The most standard example is the free scalar field theory with z = 2, known as the Lifshitz model. Anisotropic 


fixed points in interacting field theories in general can have z Æ 2 [5] [6]. Even though the original Lifshitz fixed 
points were found in anisotropic magnets where three critical lines meet, in this paper we simply define Lifshitz-like 


fixed points as any fixed points which have anisotropic scale invariance. 


of quantum criticality in condensed matter physics (see the references in [8]). See [9] {10} [TT] [12} 
for further progress on holographic aspects of this topic}? 

The other cases where 1 < p < d—1 are not only generalizations of p = 0 case but can 
also be interpreted as space-like anisotropic fixed points (see also [I3]) as is clear from the 
expression (1.4). Lifshitz fixed points with space-like anisotropic scale invariance appear in 
realistic magnets such as MnP and the axial next-nearest-neighbor Ising model [5]. They are 
also realized in models of directed percolation [7]. 

To understand holographic duals of such gravity backgrounds, it is the best to embed them 
into string theory, where microscopic interpretations are often possible by using D-branes. How- 
ever, so far there has been no known embedding of in string theory. Motivated by this 
circumstance, in this paper we will construct such anisotropically scaling solutions in type IIB su- 
pergravity. We mainly focus on the backgrounds generated by intersections of D3 and D7 branes. 
They correspond to the choice p = 2 and d = 3 and are expected to be non-supersymmetric. 
This restriction is imposed not only for the tractability of the supergravity analysis, but is also 
due to another motivation, namely to construct back-reacted D3-D7 solutions that are dual to 
the pure Chern-Simons gauge theory in the second setup of [22]. In the end, we find a class 
of solutions with the exponent z = 3/2. We also extend them to black brane solutions dual to 
finite temperature theories. 

Furthermore, we show that there exist solutions which interpolate between our anisotropic 
solutions and the familiar AdSs x Xs solutions. We also construct their numerical solutions. 
The holography suggests that our Lifshitz-like fixed points can be obtained from various four- 
dimensional CFTs including V = 4 super Yang-Mills via RG flowg®|. These flows are triggered by 
the relevant and anisotropic perturbation which gives a non-zero 0 parameter (i.e. the coefficient 
in front of the topological Yang-Mills coupling F ^ F) that depends linearly on one of the three 
spatial coordinates i.e. 0 x x3. Notice that when x3 is compactified, the perturbation induces 
the Chern-Simons coupling f A A F + 2 As as in [22], which becomes relevant in the IR. 

This paper is organized as follows. In section [| we present solutions dual to a class of 
Lifshitz-like fixed points based on D3-D7 systems with their black brane generalizations. In 
section [3] we show there exist interpolating solutions which approach the Lifshitz-like scaling 
solutions in the IR and the standard AdS; solutions in the UV. In section [4] we holographically 
calculate the shear and bulk viscosity. In section [5] we compute their holographic entanglement 
entropies and discuss how the scaling behaviors of the entanglement entropies depend on the 
direction along which the sub-systems are delineated. In section [6] we study the perturbations 
around these backgrounds and discuss the instabilities. In section [7] we present anisotropic 
solutions based on D4-D6 systems. In section [8] we summarize our conclusions. 


5 Gravity duals of another types of fixed points with non-relativistic scaling symmetry [15] have also been 


studied especially for systems with non-relativistic conformal invariance i]. 
°Here the Lifshitz-like fixed points are realized in the IR limit. It is also intriguing to consider opposite RG 


flows, where IR fixed points become relativistic and isotropic z = 1 as in [8] [0]. 


2 Holographic Duals of Lifshitz-like Fixed Points in Type IIB 
String 


In this section we will present the main result of this paper. We will construct new solutions 
in type IIB supergravity with RR 5-form and 1-form fluxes whose Einstein metrics enjoy a 
nice scaling property. Since their scaling is anisotropic as opposed to the well-known AdS5 
background, we argue that they are dual to Lifshitz-like fixed points described by certain D3-D7 


systems. 


2.1 Type IIB Supergravity 


The IIB supergravity action Sg = ss J £ in the string frame is defined by the Lagrangian 
10 
(we follow the convention in [23]) 


-2¢ 1 1- : 
L = V Zge™?® (R + 40 60™ $) — Hh; AsH- SF AFi — SPA «Ps 


1 ~ ~ 1 
ib ^P — Ca ^ Ha A Fs, (2.1) 


where F, = dy, F3 = F; — x H3, and Fs = F; — iC A H3 + 5 Be A^ F3. We set a’ = 1 therefore 
2ko = 27)". 
The fluxes obey the equations of motion: 
d* F, = *F3 A Hs, d* Fz = -H3 A^ Fs, dx Fs = H3 A^ Fs, 
d(e~?? x Ha) = Fy A xP + Fs A Fs. (2.2) 


plus the Bianchi identities: 


dH;=0, dF, =0, dFk3=H3AF,, dFs = H3A Fs, (2.3) 
and the self-dual constraint for Fs: 
« Fy = Fs. (2.4) 
The dilaton equation of motion is 
RAAV UV" oH Hun pH? —4AV yoV" o =0. (2.5) 


And the Einstein equation becomes 


1 
Run + 2VuMVno+ qouna 


24 F = ABCD 
e? Fy apcpF y ; 


(2.6) 


1 1 looks 1 
= zHmasHy’” + 50 Fu Fw + qe FMB hy AY + TEET 


where 


1 7 : 
=e Fapcpp F POPE. (2.7) 


bee o 
A = e’ ðmxð x + ao eer T 


2.2 D3-D7 Ansatz 


We start with the (string frame) metric ansatz that preserves the three-dimensional Lorentz 
symmetry SO(2, 1): 


de = err (gr a dx? 4+ dy?) ie e2h(r)+2a(r) day? pa e2c(r)—2a(r) qr? $ a A (2.8) 


We require the five-dimensional compact manifold X; to be a unit-radius Einstein manifold with 


the same Ricci curvature as the unit-radius 9%, i.e. it satisfies 
Rag = Aga. (2.9) 


The simplest example of Xs is obviously the unit radius sphere S°. The self-dual 5-form and 


1-from fluxes are given in terms of constants a and 8 by 


F; = a (Qx; + *Ox;), (2.10) 
F, = dy = pdw, (2.11) 


where y is the axion field (i.e. the RR 0-form potential) and Qx, is the volume form of X5. The 
fluxes and satisfy the equations of motion (2.2). We also assume that the dilaton ¢ 
only depends on r and both 3-form fluxes (H3 and F3) vanish. Our ansatz, which looks rather 
different from [8], is motivated in part by an attempt to construct back-reacted solutions of the 
D3-D7 intersecting systems introduced in [22], as will be explained in detail later. 

Under this ansatz, the equations of motion for the metric and the dilaton ((2.6) and (2.5)) 


are summarized as follows: 


2 
jaa _ BŽ —2a-h+3b46c,.5 | QÊ —4e+3b+h 8 
b'e] =e r +T r», 


B? a? 
(a + hen} — — Z e72a—h+3b+6c,,5 + E 
eran 4 99 Dm. h4+3bt6c OC 4c+3b+h,,—5 
(c+ logr)'e?*] = se z— ape a—h+ +6C,, rid c+3b+ rT , 
20 2 2 
(2z+c- aye") _ Ge” _ Egat _ S ne 
1 5 3 10e? 5 
27! + d! _ a ai 2(z')? a 5 (hy a3 ah! + 2(c')? + C A: a’)d $ O = = a =0. 
Here we have defined 5 5 i 
z= pot 5lgr+a+2c+5h- 4. (2.12) 


The derivative of a function f with respect to r is denoted by f'(r). An observation, which will 


be useful in the next section, is that a linear combination of the first four equations gives 


(2b — 2a — @ — 2h)'e”] = 0. (2.13) 


2.3 D3-D7 Scaling Solutions: Holographic Duals of Lifshitz-like Fixed Points 


Since we are looking for scaling solutions (namely solutions invariant under scale transforma- 
tions), we require all metric components in to be power functions of r. In other words, 
the functions a,b,c, z and ¢ are all logarithmic functions of r. For such a scaling ansatz, the 
equations of motion and reduce to algebraic equations and the solution is easily found 
to be: 


a(r) = sles — logés, b(r)= m logr +bo, c(r)= (-1 +$) logr + co, 
h(r) = ŽŽ logr +logé. + ho, 6(r) = Š logr + go, 


g= eco $0 B=A Ze eho—co— po, (2.14) 


where bo, co, ho, Qo and €, are arbitrary constants. Es corresponds to the degrees of freedom 
of the reparameterization of r, while bọ and ho correspond to the rescaling of the (t,x, y, w) 
directions. 


Without loss of generality, we choose 


1 11 11 
é= l; bo = co + 5 log Tp ho = co + log 75> (2.15) 
and the solution in the string frame reads: 
S 2 2 dr? ee Be: 
ds? = R? |r3(—dt? + dz? + dy? )+r3dw? +— | + A,r dsy,, (2.16) 
r3 

where R2 = ER? = e7°, And the dilaton scales with r as 

e? = 736%, (2.17) 


v22 


where e? = E 

Since the dilaton depends on r non-trivially, it is helpful to discuss the metric in the Einstein 
frame. Indeed, later we will see explicitly that a large class of scalar fluctuations around this 
solution can be described by Klein-Gordon equations on curved spacetimes based on the Einstein 


frame metric instead of on the string frame metric. The above solution in the Einstein frame is 


dr 
ds% = R? i 2(—dt? + dr? + dy? ) +r3dw? + l + Rds, (2.18) 
where the radii i 
P= R ge an (2.19) 
11 2 


The metric (2.18) is invariant under the scaling 


(t,£, y, w,r) > (àt, AT, Ay, AF w, =) : (2.20) 
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and therefore is expected to be holographically dual to Lifshitz-like fixed points with space-like 
anisotropic scale invariance. Note that the metric is equivalent to (1.4) with z = 3/2, 
p=?2 and d=3. 

By redefining the radius coordinate p = r3 and rescaling (t,x, y,w) accordingly, we can 
rewrite the metric (2.18) into another illuminating form 


dp? 


dst = R? |o (=d? + da? + dy?) + pèdu? +| + Redsk,, (2.21) 


This can be regarded as gravity duals of Lifshitz-like fixed points with z = 3/2. It coincides 
with the metric (1.2) with p = 2 and d = 3. 


2.4 Holographic Interpretation in terms of D3-D7 System 


Since our solution is sourced by the RR 5-form and 1-form flux (2.11), we expect 
it to be interpreted as a D3-D7 system in string theory. When we compactify the w direction 
such that w ~ w + and place N D3-brane along the (t,2,y,w) directions and k D7-branes 
along the (t,x, y, X5) directions: 


Max StxX5 |t zı y rlwis, so 83 s4 85 


N D3 |x x x Ix 
kD7 |x x x x xXx xX xX xX 


these N D3 and k D7 branes can source the desired RR 5-form and 1-form fluxes with 
(2r)4N k 
TEA = =, 2.22 
EES ma eee) 
This brane configuration is the same as the one constructed to model the fractional quantum 
Hall effect in [22]. 


The number of the D3-branes determines the radii R and R in the scaling solution (2.18): 


12 = Tt 
RP = —R? = 2,/—___N. 2.23 
11 Vol(X5) Ca) 


For Xs = S5 (whose volume is 7°), R? = 2R? = 2/aN. The number of the D7-branes gives 
the string coupling at r = 1: 

do V 22 L 

e” = —__., 

3 k 

Now we say a few words about the field theory living on this D3-D7 system. We take 


(2.24) 


X; = S° to simplify the arguments. If we start with M D3-branes, whose low energy theory is 
the four-dimensional M = 4, SU(N) super Yang-Mills theory, then the additional k D7-branes 
will source a non-trivial axion field y = tw, which in turn induces a w-dependent 0 term (i.e. 
the topological term) of the Yang-Mills theory 


T xlw)TrF A F. (2.25) 


For finite 8 and k, w-direction is compactified. After integrating over w, the 4D topological 
term (2.25) becomes a 3D Chern-Simons term at level k: 


k 


2 
=a Tr | AAF +2A?]. 2.26 
AT R12 | t 3 | ( ) 


Now we have two choices of the boundary condition for the w-circle: periodic or anti-periodic. 
If we impose the anti-periodic one, all fermions will become massive. This breaks all supersym- 
metries and gives masses to scalar fields through quantum corrections. In the IR limit, only 
a pure Yang-Mills term is left of the original 4D M = 4 super Yang-Mills part of the action. 
Since in the IR limit, the Chern-Simons term dominates this Yang-Mills term, the final three- 
dimensional theory is a pure Chern-Simons theory. In [22], this D3-D7 system was constructed 
to holographically model the FQHE precisely because it flows to the pure Chern-Simons gauge 
theory in the IR. In this model, the AdS/CFT correspondence in the IR limit manifests itself 
as the level-rank duality of the pure Chern-Simons gauge theory. 

On the other hand, if we take k — oo (and simultaneously L — oo) while keeping £ finite, the 
w-direction is non-compact and the field theory is four-dimensional. Even though the interaction 
looks non-local at the first sight, its contribution to the equations of motion is actually 
local. This remarkable property occurs only when x(w) is a linear function of w (as is the case 
here). 

One might still doubt any relations of our new background to the N = 4 super Yang- 
Mills theory as it is not asymptotically AdS5. One might also worry that the dilaton blows 
up near the boundary r — co. However, as we will show in the next section, we can in fact 
construct solutions which interpolate between the AdS; and our scaling solution (2.18). This 
interpolating solution can be considered as the dual of the RG flow between the two systems. 
Notice that this caps off the strongly coupled region of the scaling solution. We will also present 


anisotropic solutions for analogous D4-D6 systems in section [7] 


2.5 Black Brane Solutions and Entropy 


One more interesting fact about our scaling solutions is that we can straightforwardly generalize 


them to black brane solutions which have regular event horizons. The metric in the Einstein 


frame is dr? 
ds} = R? |r?(—F(r)dt? + dx? + dy?) + r3 dw? + F + Ridsx,, ele 
re? F (r 
where 
F(r) =1-&. (2.28) 
r3 


The constant u represents the mass parameter of the black brane. The dilaton and RR fields 


remain the same. 


Requiring the smoothness of the Euclidean geometry of (2.27) gives the Hawking temperature 


TH = eT. (2.29) 
The Bekenstein-Hawing entropy is then 
T? 2 mè 
Sou =7: (xs) - N? -TÈ Va- L, (2.30) 
where y is a numerical factor 
y= 23-36-1178 - 3 ~ 3.729 (2.31) 


and Və represents the area in the (x, y) direction. The entropy is proportional to N? and 
thus is consistent with the planar limit of a certain gauge theory. 

Notice that the power 8/3 of temperature in can also be obtained from a simple 
dimension counting. From the metric (2.27), the coordinate w has the fractional dimension 2/3, 


while each of (t, x,y) carries the unit dimension. 


3 RG Flow in AdS;/CFT, and Scaling Solution 


In the previous section, we find a new scaling solution of the D3-D7 system in type IIB su- 
pergravity. To clarify its physical interpretation, we will show below that we can construct 
interpolating solutions that approach the AdS; x X5 solutions in the r — oo limit (i.e. UV 
limit of the holographic duals) and the scaling solutions in the opposite limit r — 0. Then 
via the AdS/CFT correspondence, we can argue that the system dual to our scaling solution is 
connected to the one dual to the AdS; through the RG flow. 


3.1 Further Reduction of Equations of Motion 


To find the interpolating solution, we start with the general form (2.8). To simplify the problem 
we impose some extra constraints which are consistent with both the AdS; and the scaling 
solutions. 


First, we can make the function a(r) vanish by a reparametrization of r: 
alr) = 0. (3.1) 


Secondly, recall that we showed [(2b — 2a — ¢ — 2h)'e?]' = 0 for generic solutions. In fact, both 


the AdS; and the scaling solution satisfy a much stronger condition 


(2b — 2a — $ — 2h)’ = 0. (3.2) 


Since we are looking for a solution that interpolates between the AdS; and the scaling solution, 


it is reasonable to impose (3.2) as a simplifying ansatz, namely 
1 = 
h(r) = b(r) — a(r) — 59(r) + ho, (3.3) 


where hg is a constant. Similarly, since both the AdS; and the scaling solution have the nice 
property that in the Einstein frame the radius of S° is a constant, we will also impose this 


condition on our interpolating solution, namely we require 
olr) = 4c(r) + 4logr + ĝo, (3.4) 


where ġo is a constant. These relations (8.1), (3.3) and (8.4) are the constraints mentioned and 
are assumed throughout this section. 


Under this ansatz, the Einstein frame metric becomes 


2b(r)—2e(r)— $0 2b(r)—6c(r) _ ġo dr? 


zoj a 
E (-dt?+dz?+dy?) +e- 340+2h0— dw? +e Fy te dsk, (3.5) 


ds, = 
E r2 


and the equations of motion are greatly simplified: 


a= 4e?0, (3.6) 
2 23 16c 

b" = 5 22y 102 = 2E H 24b'd =e", (3.7) 
r r r 
4 14 5c 

J 24 — Ua 6y* d, (3.8) 
r r r 


1 ae ed 
7 eae ak = 6 — 8 — 6r7b? — 16rd — 8r°d? + 18rb'(1+ rc). (3.9) 
We can confirm that the derivative of the r.h.s of (3.9) is vanishing if (8.7), (8.8) and (8.9) are 
satisfied. This means that the constraint (8.9) is consistent with (8.7) and (8.8). 


3.2 Interpolating Solution between AdS; and D3-D7 Scaling Solution 


Now the problem amounts to solving the system of two coupled first-ordered nonlinear ODEs 
and under the constraint (3.9). First, notice that and involve only the 
derivatives of b and c, thus once we find a solution to them, we can simply choose the integration 
constants of b(r) and c(r) such that they satisfy the constraint — as long as it allows the 
r.h.s of to be ponove Therefore essentially we only need to solve and (8.8). 

Next we redefine the radial coordinate r and the derivatives of the functions b(r) and c(r) 


as follows 


s=logr, B(s)= (3.10) 


THowever as we will show later this requirement is actually automatically satisfied by the interpolating solution 


we are looking for; thus it does not impose any additional constraint. 
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Then the equations (8.7) and (8.8) are simply a pair of first-ordered non-linear ODEs: 


B=24+424B — 16C — 10B? + 24BC — 8C?, (3.11) 
Ò = 4+ 14B — 4C — 6B? + 14BC — 20°. (3.12) 


where B = EA A physical solution also needs to satisfy 


(9B — 8C — 8)? — (33B? + 48) < 0, (3.13) 


due to the constraint (3.9). 
The dynamical system (8.11) and (8.12) has four fixed points (B,C)* which can be classified 


into two pairs: 


7 1 
B,C)* = (+1,-1 and B,C) = ——,łŁ—— - |1}. 3.14 
(B,0}" = (#1,-1) (Bc) = (+-+) Gay 
Inside each pair, the two fixed points are related by a coordinate redefinition r — 1 thus are 
equivalent. The fixed point 
(B,C) ras, = (1, —1) (3.15) 


corresponds the standard AdS; x X5 sölutton 


Our scaling D3-D7 solution (at zero temperature) corresponds to the fixed point 


Co. = (= =a — i) ~ (1.2185, 0.8259). (3.16) 
One can easily see that the metric (3.5) with is equivalent to via the redefinition 
of radial coordinate r > r¥33/6, Since the two fixed points with “—” sign are equivalent to the 
two with “+” sign and are disconnected from them, we will not consider those any further. 
Now let’s study the behavior of this dynamical system and (3.12). Near the AdSs 
fixed point (8.15), the eigensystem of the linear perturbations (defined by B(s) = Bigg, + €(s) 
and C(s) = Chas, + €c(s) for €b, €c < 1) is 


E€ = —20e4 + 24€,, Ee = —12e, + 14e,, (3.17) 


and both eigenvalues are negative: Aj = —4 and Ay = —2; therefore the AdSs solution is a stable 
fixed point as the system flows to the UV (i.e. r — ov). 

On the other hand, the eigensystem of the linear perturbation near the D3-D7 scaling fixed 
point (defined by B(s) = Bžearing + Mo(s) and C(s) = Ccating + Ne(S) for m, Ne < 1) is 


. 116 152 . 70 94 
Nb = 733 + rag "Ie Neo == ag I T ag e (3.18) 
8The fixed point (B,C)* = (—1,—1) can be considered as its conjugate since they are connected under a 


coordinate redefinition r —> 2. 
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In contrast to the stable AdS; fixed point, this scaling solution fixed point has one negative 
— 88+ /105 v105) and one positive one (Aj = H% VI Vy; Therefore the fixed point 


corresponding to the scaling solution is unstable. Near the ie ane of this scaling fixed 


eigenvalue (Ay = 


point, there exist two special trajectories: one corresponding to the negative eigenvalue 1 and 
one to the positive àz. The fixed point behaves like a UV (resp. IR) fixed point when approached 
along the trajectory corresponding to the negative (resp. positive) eigenvalue. When the fixed 
point is approached along a generic direction, the trajectory only passes near it and then turns 
to flow to infinity — only one fine-tuned trajectory can reach the fixed point. 

Figure [I] shows the global behavior of this dynamical system and (8.12). Figure [1 (a)| 
is generated by numerical computations and the salient features are schematically highlighted in 
the hand-rendered Figure {1(b)] The arrows point in the direction from the UV (s = oo) to the 
IR (s = —ov). It is easy to see from the direction fields that the AdSs5 fixed point = os dot 
at (1, —1)) is a stable UV fixed point while the scaling fixed point (the blue dot at (—& JE 7 —1)) 
is unstable. Recall that a physical solution needs to satisfy (8.13). The allowed region in the 
flow diagram are between two hyperbolic lines given by 9B —8C — 8 = +\/33B? + 48. The black 
curve in Figure [1(a)] is the one with the “+” sign; and the other one with the “—” sign is its 
mirror image in the upper-left corner but is out of the range of Figure[1(a)| It is clear that both 


the AdS; fixed point and the scaling solution fixed point are in the allowed region. And since 
the AdS; is a stable UV fixed point and there is no critical surface separating it from the scaling 
solution fixed point, there exists a trajectory emanating from the AdS; fixed point and flowing 
to the scaling solution fixed point. Namely there exists a solution that interpolates between the 
AdS; x X; solution in the UV (r — oo) and the D3-D7 scaling solution in the IR (r = 0). 
Now to solve the interpolating solution, we first choose the integration constant. First, do 
is determined only by the 5-form flux: do = — log $ = —log (vata ). Then without loss of 


o 
=% 


generality, we set ho and choose the boundary condition for (b,c) to be 


b(r) > logr, c(r) + —logr at r — +00, (3.19) 


v33 v33 
4 . 
where bo = $ log lea) — b log 3 — $ log 3. Then the dynamical system (3.11) and (3.12) 


plus the constraint 


7 1 
b(r) = == logr +bo, cfr) > (= - i) logr+bo at r—0, (3.20) 


lB — 8C — 8)? — (33B? + 48)] = -— É eM <0, (3.21) 
64 (vis) N? 
and the boundary conditions determine an interpolating solution that approaches the 
AdSs x Xs solution 


d 2 
ds? R? |r?°(—dt? + dz? + dy? + dw?) + = + R?ds¥, 


e? = e”, (3.22) 
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(a) (b) 


Figure 1: The (B,C) flow diagram. The horizontal and vertical axes are B = b and C = ċ, 
respectively. The arrows point in the direction from the UV (r = co) to the IR (r = 0). In the 
left figure, the blue dot at E WEE] —1) ~ (1.22, —0.83) is the unstable fixed point corresponding 
to the scaling solution; the green dot at (1,—1) is the stable UV fixed point corresponding to 
the AdSs solution. The green line running through the scaling fixed point corresponds to the 
negative eigenvalue 1 while the red line corresponds to the positive one Ag. The black curve is 
given by 9B — 8C — 8 = /33B? + 48 and corresponds to a pure D3 solution; the allowed D3-D7 


solutions are above this curve. 
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in the UV and the D3-D7 scaling solution 


12 dr? 2. 
dst, = R?[rv33 (—dt? + dx? + dy”) + 2 + pr v75 dw + R?dsx,, (3.23) 
r 
e = eto etor TB, (3.24) 
in the IR. Here R? = e~90/? = 2, lychee and p? = e7% R?. And the fluxes are 
Fs = 4R* (Qx, + *Qx,), F, = pdw, 


throughout the system. 
We can easily solve (b, c) numerically for arbitrary fluxes. Systems with different flux numbers 


differ only in their speeds in approaching the fixed points. Figure [2] shows the system with 


TEE ET = 1 as an example. Figure graphs the behavior of (b, ċ), and more directly, 
Vol(xs) 


Figure [2(b)| presents the scalings of the (t, x, y)-directions, w-direction, and e? in the Einstein 
frame as the interpolating solution flows from the D3-D7 scaling solution in the IR (r = 0) to 
the AdSs in the UV (r — ov). 


=p 
0 
1- 


B,C 


WwW 


e {phi} 


3 
log(r) 


al 


(a) (B,C) of the interpolating solution. (b) Scalings of the interpolating solution in the 


log(r) 


Einstein Frame. 


Figure 2: An interpolating solution with fluxes satisfying —_——; = 1. In the left figure, 
‘Volo 


1 — 1) ~ (1.22, —0.83) in the IR to (1,—1) in the UV. The right 


(B,C) = (6, é) flow from (Ja Ka 
Einstein 


figure shows that, in the frame, the scalings of the (t, x, y)-directions, the w-direction, 


and e? flow from (=, =, — 1.04, 0.70, 0.70) in the IR to (1,1,0) in the UV. 
TR V33? T 
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3.3 Interpolating Solutions as Holographic RG Flows 


Here we try to interpret our interpolating solutions via AdS/CFT. To make the argument 
simple we focus on X5 = S°. Consider the standard AdS;/CFT, for the M = 4 super Yang- 
Mills. We can perturb the M = 4 super Yang-Mills by many relevant operators O1, O2,- as 
65 = f dx*[giO1 +9202 +--+]. 

Our solution which flows from the AdS; x S3 to the scaling D3-D7 solution is dual to an 
anisotropic RG flow triggered by a non-supersymmetric relevant deformation (2.25), called O1. 
This RG-flow eventually ends at the IR fixed point that is dual to the scaling solution. Therefore 
O, becomes irrelevant in the IR limit. However, our sketch of the holographic RG flow (see Fig. 
tells us that this flow is unstable and signals the presence of another relevant operator O2 
which becomes relevant even at the IR fixed point. So if we slightly perturb this flow by Og, 
then the RG flow passes near the IR fixed point and eventually goes to infinity. Its asymptotic 
behavior will be derived analytically in the next subsection. However, we can still fine tune such 
that there is no Og generated (i.e. g2 = 0) to realize the IR fixed point governed by the scaling 
solution. This is exactly what experimentalists usually do to realize an unstable fixed point 

A similar situation occurs when we perturb the Heisenberg model by an anisotropic Ising-like 
interaction: H = —J lua (oroj + ojo} + 070%) —D Jiz) % 97, Where D > 0. In this analogy, 
the UV fixed point is the Heisenberg model fixed point and the IR one is the Ising model fixed 
point. This structure of the RG flow is rather generic and is one example of the phenomenon 


called crossover [7]. 


3.4 Exact Solutions without D7-brane Charge 


In our interpolating solution, the flow from the UV AdS; fixed point reaches the IR D3-D7 
scaling solution fixed point only if it starts along a specific direction. One wonders where the 
flow would end up if its initial direction slightly deviates from the fine-tuned one. A closer look 
at Figure [1(a)] tells us that if it deviates from the desired direction slightly to the left, it will 
turn back before reaching the scaling solution fixed point and asymptote to the B < 1 part of 
the black curve given by 9B — 8C —8 = 33B? + 48 (due to the constraint (8.13)). On the other 
hand, if it deviates slightly to the right, it will pass near (but not hitting) the scaling solution 
fixed point and then bend slightly downward and finally asymptote to the B > 1 part of the 
same black curve. 

The curve given by 9B — 8C — 8 = V33B? + 48 is the solution interpolating between the 
AdS; in the UV and some other solution in the IR. Since it saturates the inequality (8.21), there 
is no D7 brane charge (3 = 0): it is a solution of the pure D3-brane system. Not only is it 
important because it gives the asymptotic form of what our interpolating solution decays into 
Tt is feasible because as long as the starting point is close enough to the critical line (which corresponds to 


the fine-tuned trajectory that hits the fixed point), the trajectory will stay for a very long time near the fixed 


point to allow the measurements [7]. 
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when perturbed by a relevant operator, it is also interesting as a pure D3-brane solution other 


than AdSs5. In this subsection, we will solve it analytically. 
First, eliminating C from equations (3.11) and (8.12) gives 


11 
— =44 —R?_ 
ds t7 


which can actually be solved analytically. There 


-B 
4 


3 


48 + 33B2, (3.25) 


are two solutions distinguished by the + sign: 


e + 6 get +1 


B.(s) = e8s = 1 ? 
e85 + 44/ pets —1 
C+(s) = BL] (3.26) 
which lead to 
> 3/3, +1 1 R 
b+(s) =b -s7 n® As i+q (e° — 1), 
. 1/3, e#+1 
c+(s) =O) — Sse 3 Tl log Asc] (3.27) 
We can choose ho = $o and bọ = Č = 0; then the metric in the Einstein frame and the dilaton 
are simply 
ž o| /r8—1 /rf+1 t3y i : 4 P r8—1 /rf+1 Taya 9 
dh = R “a (==) (—dt* + dx* + dyf) + a (4) dw 
d 2 
+ 2 + ds, : 
3 
4 1 +2 TT 
e? = (5 + 7) efo, (3.28) 
5 2 pz 


Both solutions in become AdSs x Xs in the UV fixed point (r — oo); however when 
going towards the IR, they become singular at r = 1. The solution with “+” (resp. “—”) sign 
covers B > 1 (resp. B < 1) part of the curve 9B — 8C — 8 = V33B? + 48. Even though these 
zero temperature solutions become singular ar r = 1, their corresponding black brane solutions 
at finite temperature are expected to be smooth. An explicit construction of such black brane 
solutions is left as a future problem. Also it would be very interesting to understand what 


anisotropic relevant deformations of the NV = 4 Yang-Mills are dual to this background. 


4 Hydrodynamics 


As we obtained the black brane solutions (2.27), we can consider the hydrodynamic behavior of 
their dual field theories from the supergravity side [24]. In this section, we especially determine 
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the shear and bulk viscosities from dispersion relations for the corresponding quasinormal modes 
of fluctuations around the background (2.27), following [26]. 

For this purpose, we first reduce S5 part and consider fluctuations around the resulting 
five-dimensional background. The procedure of the reduction and derivation of the linearized 
equations of motion for the fluctuations are summarized in Appendix[A.8.1] In these equations, 
fluctuation of the resulting five-dimensional metric Sg = H, dilaton 6¢ = ọ, axion dy = n, 
the trace part of the S° metric 7 and the five-form flux dF; = fs appear generally. 

In this section, we assume that the w direction is compactified and the neglect the momentum 
in this direction so that we extract the effective 2+ 1 dimensional holographic dual theories from 
the total 3+ 1 dimensional ones. In this situation we can choose the momentum in y direction. 
By considering the symmetry of the background, we can decompose the fluctuation into channels 


decoupled from each other in the linear order: 


shear channel Hiz, Hry, Her, (4.1) 
scalar channel Hzw, (4.2) 
sound channel Hu, Hiy, Hir, Hex, Hyy, Hyr, Hww, Hrr, 9, T, ftzywr, (4.3) 
“axion channel” Heiss Hyw, Hrw, 1- (4.4) 


In this paper, we consider the shear and sound channel only, since we can read off the shear and 
bulk viscosity from dispersion relations for the quasinormal modes of these channels. 
Before deriving the dispersion relation, we need to consider the gauge fixing of the fluctua- 


tions. Let us assume that the fluctuations are of the form 


Hw = A (rje em), p= (re ieta), 
n = n(rje En), mE nhre), (4.5) 


Fuvpor = Fuvpord (re tea) : 


Especially, for later convenience, we rewrite some of their components as 


Hy(r) E R?2Ha(r), Ha(r) = ReHalr), (4.6) 
alr) R? Hig(r), Huw(r) = R2c4/? Huy (r), (4.7) 


where a,b = x or y, and c? = r?F(r) and c2 = r?; the function F(r) is defined by (2.28). 
As in Appendix we then take the radial gauge H,, = 0 in the following discussion. 
Here we notice that there still exist residual gauge degrees of freedom under an infinitesimal 


diffeomorphism z” > x! + ¿r with €¢ = €#(r)e ttm); 


Hy f Hw Vis Vu ea ce E O" Q, (4.8) 
N => N= EOX, T >T Ep" Jo - (4.9) 
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Here the covariant derivative is defined by the background metric. In [25], it has been noticed 
that one can derive the dispersion relation for the quasinormal modes by defining gauge invariant 
quantities, instead of performing the gauge fixing completely. Following this, we define the gauge 


invariant quantities as 


Zı(r) = qHiz(r) F wHgy(r), (4.10) 
for the shear channel and 
oC 2 3 (9 CeCe 2 

Zo(r) =4 z Hulr) T 2qw Hy, (r) +w Hyy(r) + 5 q ae w“ ) (Aee(r) + Hww(r)), 

(4.11) 
1 

Zo(r) = g(r) — 5Hww(r), (4.12) 
Z,(r) = 1(r), (4.13) 


for the sound channel. Below, we derive equations for these quantities by using linearized equa- 
tions of motion summarized in Appendix|B]and then the dispersion relations for the quasinormal 


modes by imposing an appropriate boundary condition. 


4.1 Shear Viscosity 


From the equations (B.1),(B.2) and (B.3), we obtain the equation for Z, as 


2 2 2.7 of 2 
24 SFR eT Go a 
where the dimensionless frequency w and momentum q are defined by 
= ae q= re (4.15) 
respectively. 
In order to solve this equation, we consider the asymptotic behavior of Z1 first. As r — oo, 
Z, ~ Ar? + Br~s with constants A, B generally. For the quasinormal modes, A is set to 


zero, or, in other words, we should impose the Dirichlet boundary condition Z; = 0 at r = oo. 


As r > uI, Zı ~ (F(r))*™/2 and, in order to ensure that only incoming waves exist at the 
horizon, we take the one with negative sign. 


By taking the hydrodynamic limit q « 1 and w « 1, we can perturbatively determine Z; as 
. j 2 
Z= C(F (ry) (: + oF (r) + O(q?, w) ) ; (4.16) 


where C is a constant. Here we assumed that q? and w are of the same order as usual for the 
shear channel. Then by imposing the Dirichlet condition 7, = 0 at r = œœ, we find w = a, 


Comparing with the hydrodynamic relation for the shear channel 


we can find that the universal bound for the shear viscosity 7 to the entropy density s ratio [27] 


is saturated for the current case: 


9 (4.18) 


4.2 Bulk Viscosity 


From the equations derived in Appendix [B.2| we obtain the differential equations for Z,, Zr 
and Zo as 


1 w2 a 
Zy + In e A + 2 (= - 5 ) Z 


t t Cy 
44 11(u3 — 2u2r11/3 + 2r11 
W E i L E (4.19) 
9c? 2728/3 (711/3 — u) 
1 [w af 32 
Z! +W (EEZ | Ss), = Zr = O, 4.20 
T n (ci Cy ) T + e 2 c2 c ( ) 
Zo +F(r)Z+G(r)Zo + H(r)Zy = 0, (4.21) 


where 
625p°q* — 3584(q? — w?)?rlt 
220p (1574? — 112w?)r 1/3 + 32u(149q? — 24w?) (q? — w?)r??/3 
3r(u— re) au = 10r a) l 

3025u4q* = 5u?q?r 16/3 (3334? — 1936r?) (q? — w?) = 2304r38/3 (q? — w?)3 

O gru- rB ue 6r ABe- 
3744ur? g (q? — w?)? + 5p3q2r?/3(45q4 — 2541q2r? + 1936r?w?) 

22(1536r1! (q? — w?)3 — 1440ur??3 g?(q? — w?)?) 

22(— u? (4254? — 550w?)q* + 20y2q2r!!/3(83qg4 — 171q2w? + 88w*)) 
AB? (ue — riB) 5u — 16r4/3(q@2—w2))2 


Then, as in the case of the shear channel, all we have to do next is to solve these equations 


F(r) 
(4.22) 
G(r) 
(4.23) 
H(r) = — 


. (4.24) 


by imposing the incoming wave boundary condition at the horizon and the Dirichlet boundary 
condition at infinity to derive the dispersion relation for the quasinormal mode of the sound 
channel. We also take the hydrodynamic limit and assume that q and w are of the same order, 
as is expected for the sound channel. For Z, = 0, we can see that the nonsingular solution for 
is a constant, which turns out to be zero due to the Dirichlet boundary condition Zy = 0 
at the boundary r = oo. By substituting this into (4.19), we obtain Z = 0 in a similar manner. 
Then from these results and (421), by using a similar ansatz Zo = (F(r))~*"/*Yo(r), we can 
perturbatively determine Zp as 


Zo = Č(F (r) 2 (1 sO Pu ren) + O(@?, 0°, aw)) | (4.25) 


where Č is a normalization constant. As a result of the Dirichlet boundary condition Zo = 0 at 


r = œ, we obtain the dispersion relation for the sound channel 


iB os 


Let us recall a hydrodynamic relation for the sound channel in the noncompact (d+1)-dimensional 


spacetime 
= n 2 
v=o- izh (+E) +, (4.27) 
H 


where c, and Ç are the sound velocity and the bulk viscosity, respectively. Since there are 
two noncompact spatial dimensions for the dual field theory now, d = 2 for the current case. 
Therefore, by comparing (4.26) with (4.27) with d = 2, we obtain 


3 
= 7 2= (4.28) 


As for ¢/n and c, of strongly coupled gauge theory plasma in d noncompact spatial dimensions, 


it is conjectured in that they satisfy an inequality 


Ç l 2 
7 = 2 G = 2) : (4.29) 


In our case, this inequality is saturated. 


5 Entanglement Entropy of D3-D7 Scaling Solutions 


When a quantum system is divided into two subsystems: A and its complement B, the von 
Neumann entropy S4 = —Trp, log pa (where py is the reduced density matrix after tracing out 
B) is called the entanglement entropy. The scaling behaviors and certain universal_4 coefficients 
of the entanglement entropy encode important information on the degrees of freedom and non- 
local correlations of the system [29] [30]. 

For an anisotropic system, an interesting question is “how does the scaling behavior of the 
entanglement entropy depend on the direction along which the subsystems are delineated?” In 
this section, we will study the entanglement entropy of various subsystems A of the (x,y, w) 
space at the boundary (r — oo) of the 5D part of the D3-D7 scaling solution (2.18). The field 
theoretical computation of the entanglement entropy is expected to be difficult as the system 
will be strongly coupled. We will instead compute its holographic dual on the gravity side. The 
holographic dual of the entanglement entropy of a subsystem A is given by 


Area™” 
T 5.1 
SEE TE (5.1) 


‘Here “universal” means the independence from the choice of different delineations of subsystems. 
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where Area” is the area of the three-dimensional minimal surface that lives inside the (r, x, y, w) 
space and borders on the boundary OA of the subsystem A [BI]. 
After a coordinate transformation r = L, the metric of the D3-D7 scaling solution (2.18) 


becomes 


dst = R2 (= + dx? a. + dz? 5 = 
ž 


) + R?ds.. (5.2) 
For X; = S°, R? =2V/7N and R? = H R?. We consider the full boundary system given by 
€ [0, Le], y € [0, £], we [0, Lal. (5.3) 


Among the various types of subsystems, we will only study the easiest types: stripes with either 


x or (inequivalently) w restricted to a smaller length. 


5.1 Entanglement Entropy for Subsystem along z-direction 
Let’s first consider the subsystem A cut out along the x-direction: 
€ [0,2 < Lz], y € [0, Ly], w € [0, Ly]. (5.4) 


The three-dimensional minimal surface bordering on OA is given by the embedding function 


z= z(z): 


ao se 
; Ly e a ae R z fs 1 
Area” - f" d f dy |. dw— -23 ) + (2+) Siite dxGvi+ 2%, 
0 


(5.5) 


where d=1+4+1+ 2 = g is the total scaling of the boundary system. 
This is a Lagrangian system with £ = FR3LyLyaVv1 +2. The z(x) that minimizes the 


surface area is then given by the equation of motion 


= ewer (5.6) 


zZ 


& 
I 
Q 


where z, is the peak of z on the minimal surface, across which z’ changes sign. It can be solved 


from ¢, = 2 Jo - a which gives z, = gr. The minimal surface is then 


a 
is 1 
= 3 ; £ 
Area” =2R bylu ri -In with In= f duza To Td . (5.7) 


In has a UV divergence at z — 0. Imposing the UV cutoff by choosing the lattice spacing 


a for the boundary system, we get 


ee ae [aon E | | (5.8) 
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Then plugging the value of z, and the five-dimensional Newton’s constant G; = ge with 


Gio = 87°08, we finally obtain the holographic entanglement entropy for the subsystem divided 


out along the x-direction: 


d 
é oO frye r ming! 1 (2\%% (al) 
aE NIB} Vol(Xs) ead am (E Te 


NLL | Be - ginl (5.9) 
with d= $. 71 and J2 are numerical constants. 

Now let’s interpret the result. First, the holographic entanglement entropy is proportional to 
the area of the boundary of the subsystem 0A = L,L,, — as expected from the “area law” 
for the entanglement entropy from direct field theory computations. Second, the first term of 
the holographic entanglement entropy diverges and is cutoff-dependent. The scaling of a is given 
by the total scalings of the y and w directions relative to that of the t-direction: (1+ 2) f= 3. 

The second, finite term of the holographic entanglement entropy is more interesting: it 
is cutofft-independent therefore can be compared with the field theoretic computation of the 
entanglement entropy. Its coefficient y2 gives a measure of the total degrees of freedom. The 
scaling of y is simply the total scalings of the y and w directions relative to that of the z- 
direction: (1 + 4)/1 = 3. Since the scaling of the z-direction is the same as that of the 


t-direction, the exponents of a and ły are the same. 


5.2 Entanglement Entropy for Subsystem along w-direction 


The next easiest subsystem we can consider is to divide along the w-direction: 
w E [0, lu < Lol. (5.10) 


The three-dimensional minimal surface bordering on OA is given by the embedding function 


z = z(w): ; 
. 2 w 1 
Area? = FL aLy | dw— vV 27/3 + 22. (5.11) 
0 Z 


Then we could follow the line of the previous subsystem along the x-direction. The computation 
is straightforward but more complicated so we present here instead a simpler route which utilizes 
the result for the z-direction subsystem. 


The coordinate transformation 


aire 
z=, (tz, y, w) = 5(Ë 9, 0), (5.12) 
results in the metric 
3 -\° (-dP +d +d}? dP dv 
2 2 2 
ds? = ($2) pe + FI + T) + R*dsx,. (5.13) 
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Thus we can simply use the result from the x-direction case, with d = $ replaced by dy = 4. 


First, we write down the dictionary between values in the original coordinates and the new one. 


1. In the new coordinates, the full boundary system is 


2 2 2 
TE [0, złe], y € (0, glu We (0, gLul; (5.14) 
while the subsystem A is 
2 2 2 2 
FE [052s], JEL], WE [0, Zlu < FL. (5.15) 


2. The lattice spacing in the new coordinates is related to that in the old coordinates by 


wl 


a=a 


(5.16) 


3. The turning point of z is Z = 2 Zes 


Using the result from the x-direction case, we find the minimal surface in the new coordinates 
to be 


(5.17) 


Area” = 3R? La Ly ; 1 | 1 VT Erm 5) 
w 


After translated back into the original coordinates, it gives the entanglement entropy of subsys- 


tem along the w-direction 


dw 
P Eo 2 eee 2 OE ae alee as) 
Bee NAG Vol(Xs) YDI |e? (2 lo ri) 


Zdu 
/ 
V1 V2 
= N?’L,L — 5.18 
where D-1 = Aet = 2. The negative finite part has the same form as the result for the 


subsystem divided along the x-direction with d = 8 replaced by dy = 4 

Now let’s compare this result with the one for the subsystem along the x-direction. The 
essential features are the same. It is proportional to the boundary area L,L,. There are one 
cut-off dependent, divergent term and one cut-off independent, finite term. The scaling of the 
cutoff a is given by the total scalings of the x and y directions relative to that of the t-direction: 
(1+ 1)/1 = 2. The scaling of the Zw is given by the total scalings of the x and y directions 
relative to that of the w-direction: (1 + 1)/(2/3) = 3. Unlike the case for the subsystem along 
the x-direction, since the scaling of the w-direction is different from that of the t-direction, the 
exponents of a and ¢y are different. 


“This might be understood as follows: when all scalings are normalized with respect to the w-direction, then 


the total scaling of the (x, y, w)-space is 
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6 Perturbative Analysis 


In order to know the details of the holographic dual field theories, a basic thing to do is to analyze 
the perturbative spectra around their supergravity solutions. This offers us the information on 
scale dimensions [2] [3]. For example, scalar perturbations are dual to scalar operators O; in 
the dual field theory. These perturbative modes in supergravity are described by Klein-Gordon 
equations with various masses in the curved spacetime. Since we have the nontrivial dilaton in 
our D3-D7 scaling solutions (2.18), it is not clear a priori whether the Klein-Gordon equation 
should be obtained from the string frame metric or the Einstein frame metric. Actually our 
scaling property of the gravity solutions is only available in the Einstein frame. Also the 
study of the perturbative spectrum is necessary to determine the stability of the background. 
Motivated by these, below we will examine the perturbations around our scaling backgrounds 
(2.18). 


6.1 Description of Perturbations 


Let us analyze the perturbations of bosonic fields around the D3-D7 scaling solutions defined 
by (2.18), (2.17), and in the Einstein frame . We will closely follow the analysis of 
AdS; x S° in [32\/33]. We will denote the total ten-dimensional coordinates by M, N,.. = 0,1, 2, +: 
-,9. The five-dimensional Lorentzian spacetime (called M5) is described by the coordinate 
L, v,- =0,1,2,3,4 and the five-dimensional compact manifold X5 by a, (,--- = 5,6, 7,8, 9. 
As is clear from the IIB supergravity action in our background, the 3-form fluxes F3 and 
H3 are decoupled from the other fields (i.e. the metric, the dilaton, F5, and F) thus we can 
concentrate on the latter ones. Then the Lagrangian in the Einstein frame is written as follows 
L=/-g (r = Tetau = Zameo ó = Ta EuNrorF™NPOR) i (6.1) 
To make analysis more tractable, let us assume X5 = S° here. Then we can define the scalar, 
vector, traceless symmetric, and antisymmetric spherical harmonics on S° by Y1, Y/, Y/ 


(a 


x Bp respectively [32]. Using these spherical harmonics, the metric perturbations dgyjn = hyn 


a)? and 


can be decomposed as follows 


I 


we. o O l a 


w) pra 


where (a8) denotes the traceless symmetric part. We also denoted all indices of the spherical 


harmonics simply by J. We fix the gauge by requiring 
V Rap) = V“hpa = 9. (6.3) 
The perturbations of the dilaton and axion are defined as follows: 
6g=y'¥!, x= Y]. (6.4) 
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Finally, the perturbation of the 5-form flux F5 = dC, can be express as follows [32]: 


— pl I — pI I — pl I 
Cpm poppe g bi wos pa i Cmuznza a Dies pops Ya?» Cha p2oras E bia u2 Ylona] 


JA I I I 
Cuara203 = bl Eai aza OV p oa Caraz0304 =b Earazazad VEY . (6.5) 


They satisfy the gauge fixing condition V“Ca... = 0, and the self-duality of F5 allows us to 
. . F I 
eliminate b, popsua 800 bja uous: 
Next we substitute (6.2), (6.4), and (6.5) into equations of motion of (6.1) and derive the per- 
turbative differential equations. We omit the details of analysis here and put them in Appendix 
[A]as many parts of the calculations are essentially the same as those in [32]. 


In the end, we find that the following modes 


Scalar modes : ¢/, (Aix, bi 
pI pl 
Vector modes : (B, 6;,), 


Tensor modes : bis (6.6) 


satisfy free massive field equations which are precisely the same expressions as those in the 
AdS; x S5 badkeround|(4 In the above, the fields in the same parenthesis mix with each other. 
Therefore, these perturbations obey free field equations of motion constructed from the 
Einstein frame metric instead of the string frame metric. 

As an example, let us concentrate on the scalar mode ¢!. Its equation of motion is written 


as 


2 


where LJ, and U1, are the Laplacians of the Lorentzian part M5 and the sphere part 5°, respec- 


tively. Using the eigenvalues of Ya gy» we eventually obtai 


( SO) oto, (k = 2,3,4,- -). (6.8) 


On the other hands, the other modes %7, n/, and Ay mix with each other and obey equations 
of motion more complicated than those in the AdS5 x S° case (see Appendix [A). 


6.2 Scaling Dimensions and Stability 


We have observed that a large class of supergravity modes (6.6), though not all of them, satisfy 
the conventional free field equations with various masses via the Kaluza-Klein compactification 

12 Tn other words, for these modes the differences from AdSs x S° only come from the background metric which 
is employed to write down the free field equations. 


13 Please distinguish the total angular momentum k of the spherical harmonics of S° from the number k of 
D7-branes. 
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on $°. The scalar modes in (6.6) satisfy the equations of motion of the form (we denote such a 
scalar field by ® here) 


(O, — m°) = 0, (6.9) 


where the Laplacian Og = g” V „Vv is constructed from the Einstein frame metric (2.18). 


Now let us consider a scalar field ® on a slightly generalized scaling background 


(6.10) 


-z [d2  —dP+deidy’ dw? 
2_ p2 
ds = R (B A, sr) 


where the scaling exponent v is related to the scaling exponent z in (1.4) by v = L, Especially, 
our scaling background (2.18) corresponds to v = 2. The equation of motion is written as 


follows: 


2 2 R2 
Gl a. i + (= +p? —w? +p, 5 = 0. (6.11) 


Here w, p, and py are the frequency, the momenta in (x, y)- and w-direction, respectively. 


After redefining the wave function by ®(z) = oF Vz), we obtain the Schrodinger-like 


equation 
—W" 4 V(z)¥ = 07D, (6.12) 
where avy? 
m? Ř2 + BtvP=1 
V(z) = -r +p? +p? ge, (6.13) 


When z is small, the third term in (6.13) is small compared to the first term, assuming v > 0. 
Thus, as in the AdS/CFT case (i.e. v = 1), we can expec [4 that the stability against the 
(normalizable) perturbations is the same as that of Schrodinger problem with the potential 
-3 
V(z)= = It is well-known that the latter system is stable iff A > 0. 
In this way, we speculate that in the background (6.10), the stability condition requires 
2 

mR? > ss (6.14) 
Notice that if we set v = 1 in (6.14), we obtain m? Ř2 > —4, which is the well-known Breitenlohner- 
Freedman (BF) bound of AdS5. 


This condition can equally be implied from the behavior of the scalar field near the boundary 


z—0 
elz) ~ Az4t+ + BzA- +... (6.15) 
where 
3 3)? i 
Kee S + wey + m2R2. (6.16) 


MWe are very grateful to Andreas Karch for illuminating explanations on the stability analysis in Poincaré AdS 


spaces. 
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The quantity A+ is holographically interpreted as the scaling dimension in the dual anisotropic 
scale invariant theory. Notice that the condition (6.14) requires that the scaling dimension 


should be real-valued. 
Now let us go back to our D3-D7 scaling solutions (2.18). In this case we obtain the stability 
condition by setting v = 2/3 and R? = H R? as follows 
m? R? > s 
3 
We can apply this condition to the scalar modes in (6.6). As is clear from (6.8), all the scalar 


modes ¢! satisfy this condition. However, we find that one of the infinitely many mixed modes of 


(6.17) 


(h/,7/,b") actually has the largest tachyonic mass m? R? = —4, which saturates the BF bound 
of AdSs. This occurs only for the second spherical harmonics k = 2. Even though this tachyonic 
mode is stable in the AdS; spacetime, it seems to become an unstable mode in our D3-D7 scaling 
background as is violated 

Nevertheless, we can replace S5 with an arbitrary Einstein manifolds X5 with the same Ricci 
curvature, keeping the same scaling solution (2.18). Define the eigenvalues A of Laplacian of a 
scalar function Y such that — R? yY = AY. In this case, if 


A A 11 
ply) — pd Se 


TA i 2 — FR? (6.18) 


is always satisfied, then the above lowest mass mode becomes stable. In other words, if there is 


no eigenvalue between 
37 — 8/3 
3 
then the background can be perturbatively stable. Notice that A = 12 saturates the BF bound 


and it is the unstable mode when X; = S°. It is intriguing whether there exists such a stable 


37+ 8V3 


(~ 7.71) <A < (~ 16.95), (6.19) 


(Sasaki-)Einstein manifold. 


7 D4-D6 Scaling Solutions 


In type IIA string, the closest analogue to the previous D3-D7 system is the following D4-D6 


system: 
Max T? x Xa: |t x y rlw, wolsi s2 83 84 
N D4: |x x x x x 
k D6: |x x x X XXX 


Here (w1, w2) span a two-manifold that supports the D6 flux (we will choose it to be T? for 


simplicity) and (s1, s2, $3, $4) span a four-dimensional Einstein manifold X4 with the same Ricci 


15A perturbative instability has also been noticed in [21] for type IIB backgrounds dual to the non-relativistic 


CFT when the 3-form fluxes are vanishing. 
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curvature as that of S4. However, as we will see, this system does not support a scaling-invariant 
solutions. In this section, we will present a gravity solution that is closest to a scaling solution: 
under (t, £, y, T, w1, w2) —> (At, Ax, Ay, 5, Mwy, N3wy), the line element ds? > \~3ds?. We will 
also give its black brane generalization. 


The fluxes given by these N D4-branes and k D6-branes are: 


(27)k 


F = T2 


1 
dw, ^d Fy = (2r) N————— Q Bo =0 7.1 
WI W2, 4 ( T) Vol(X4) X435 2 ; ( ) 


where X; is a unit-radius Einstein 4-manifold (whose Ricci tensor satisfies Rij = 3gij) and Qx, 
is its volume-form; L is the periodicity of w;. This flux profile satisfies the flux equations of 


motion. The corresponding (string frame) metric ansatz is: 


ds? — eB) ( ae + dz? $ dy?) ats e 2A") Gp? +4 e2T(r) (dw? + dw?) + elds}, 


(7.2) 
with a possibly non-constant dilaton ¢(r). 
For a scaling solution, {A, B,T, Z,¢} are 
A(r) = a, logr + ao, B(r) = bı logr + bo, T = ti logr + to, 
Z = zi logr + 20, olr) = ns logr + ġo. (7.3) 


There are one equation of motion from the dilaton and five from the gravity part. For the scaling 


ansatz, they all reduce to algebraic equations and the solution is easily found to be: 


1 1 8r? 1 17 
A = (1l-—->n;)]l — log ( —————e®N | — =] — | = log ns, 
(r) D og(r) — 3 log (sa ) 5 os (J) ogn 
5 
B(r) = gs losr + bo, 
2 1 87r3 1 Qn 
T = Ín,l -log | ————e®" N -log | ek 7.4 
(r) 3 ogr + glos (Se ) + Flow (Fe i; ce 
1 1 8r? 
Z = `n;l log | ————e®N |, 
(7) = inlogt) +3 o (r ) 
e(r) = nslogr + do. 


(ns, bo, do) are three gauge parameters. 1, corresponds to the gauge freedom of r — r“, bo 
corresponds to that of rescaling the (t, x,y) directions, and ¢o gives the string coupling at r = 1 
thus corresponds to rescalings of the r and T? directions. 


Without loss of generality, we choose 


1 8r? 1, 34 1 LÊN 3, 34 
5=2, ms b N, s wle he 
ú ae ee (K ) t3 g bo = 5 log (scare) tag 


(7.5) 
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The solution in the string frame is 
2 p2 | 42 2 2 2 dr? B74 2 2 2.44.2 
ds; = Rg |r? (—dt* + dx* + dy”) + + +173 (dwz + dwz)| + Ror3dsy,, (7.6) 
r3 


8r? 


with R2 = Fe = Na H : NE. In the Einstein frame, it is 
4 
~ dr? 5 
dsh = R? [ia + dz? + dy?) + 2: + r3 (dw? + aw’) + Rr dsX,, (7.7) 
r3 


i R n? \ł/17r? L? 
with SaR = a (A ey 


This background is no longer scaling invariant. Under the scaling transformation 
r 2 2 
(t,£, Y, r, w1, w2) > (At, Az, Ay, 5p A3 wi, AS wa); (7.8) 


the metric scales as 
ds? > \~3ds?, (7.9) 


instead of staying invariant. This is not surprising since the D4-brane theory is not conformal 
in the first place. 


Generalizing to finite temperature, the corresponding black brane solution (in the Einstein 


frame) is 
2 
ds = R? r3(—F(r)dt? + dz? + dy?) + = + 73 (dw? + dw3)| + Rrsds%,, (7.10) 
F(r)r3 
with 
F(r) =1--4, (7.11) 
Ps 


where p is the mass parameter of the black hole. Its Hawking temperature is 


17 3 
The Bekenstein-Hawing entropy is 
8r? Va LA 
S=y7-(|———-]} TR- N-V, 7.13 
4 (a) H e ota 


with y = 2% -373 1773 -r3 and V is the area in the (x,y) directions. 


8 Conclusions and Discussions 


In this paper, we presented a class of gravity duals of Lifshitz-like fixed points in type IIB su- 
pergravity. They represent backgrounds with intersecting D3 and D7-branes and their Einstein 


frame metrics (2.18) enjoy anisotropic scale invariance. We also extended them to black brane 
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solutions dual to finite temperature theories. Moreover, we showed the existence of solutions 
which interpolate between our anisotropic solutions in the IR and the familiar AdSs5 x X5 solu- 
tions in the UV. Then the holography asserts that our Lifshitz-like fixed points can be obtained 
from various four-dimensional CFTs including M = 4 super Yang-Mills via RG flows. These 
flows are triggered by the relevant and anisotropic perturbation which makes the 6-angle (the 
coefficient in front of the topological Yang-Mills coupling F ^ F) linearly dependent on one of 
the space-like coordinates i.e. 6 x w. When w is compactified, the perturbation induces the 
Chern-Simons coupling f A A F + 2 AP. This theory itself seems an intriguing model worth 
pursuing in a future work, as the equation of motion becomes local in spite of the violation of 
the Lorentz invariance. 

Employing our supergravity solutions we studied the thermal entropy and the entanglement 
entropy to measure the degrees of freedom of the holographic dual theories. We found character- 
istic scaling properties in both quantities. We also holographically computed the shear and bulk 
viscosities. A more general analysis of hydrodynamics with the momentum in the w-direction 
taken into account was left as an interesting future problem. 

Moreover, we performed a perturbative analysis around our solutions and found that a large 
class of scalar modes obey the Klein-Gordon equation in the curved spacetime which has the 
expected scaling property. Also we found an unstable scalar mode when the compact manifold 
X; is S5. Since this unstable mode occurs only for a ‘d-wave’ spherical harmonics, S% might 
decay into a less symmetric Einstein manifold and be stabilized. Thus we have reason to hope 
that there exists a class of (Sasaki-)Einstein manifolds with which our scaling solutions become 
stable. Even the background with X5 = S° is still useful at least in capturing qualitative 
properties of gravity duals of Lifshitz-like fixed points, with unstable modes simply neglected. 
The construction of manifestly stable and non-dilatonic gravity duals of Lifshitz-like fixed points 
still remains as a very interesting future problem. The three-form fluxes which we assume to be 
vanishing in our solution might play an important role. 

It is also intriguing to generalize our solutions to other values of p and d in (1.2). For 
example, the simplest case p = 0 deserves particular attention. It can be formally obtained from 
our solution (2.18) by the double Wick rotation t > iw and w — it. However, the axion field 
x and its flux become imaginary-valued therefore the solution is not physical in the ordinary 
type IIB supergravity. A slightly better example which realizes the case p = 0 is a background 
based on D3-D5 systems, where D5-branes are regarded as the baryons [34]. We cannot get any 
consistent solution if we restrict to the ordinary type IIB supergravity because the tadpole for 
the H-flux is generated by F-stings which attach to D5-branes and stretch into the boundary 
[34]. To construct a solution in this background we need to add the F-string action as an extra 


term to the type IIB supergravity. Under this slightly unusual assumption, we can indeed find 
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the following black brane solution in the Einstein frame [35] 
dp 
p? A(p) 


et) = eop,  A(p) =1— Ean i? = 1027, (8.1) 
p 


ds? = —A(p)p dt? + p? (da + dy? + dz”) + L? + Ld’, 


with constant 3-form fluxes H3, F3 and the RR 5-form F;. A derivation of this solution is briefly 
reviewed in Appendix[C] At zero temperature, this corresponds to the metric with p = 0, 
d=3andz=7. 

Finally, it is also intriguing to apply our backgrounds to realistic condensed matter systems. 
Our D3-D7 model was originally introduced to model the holographic dual of fractional quantum 
Hall effects in string theory (refer to for other holographic realizations of quantum 
Hall effects). Therefore one of the future problems is to calculate physical quantities such as 
finite temperature corrections to the Hall and longitudinal conductivities in this theory. The 
anisotropic critical points we found in this paper may also be useful to analyze the systems like 


liquid crystals and some anisotropic spin systems. 
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A Perturbative Analysis 


In this appendix, we consider perturbative fluctuations around the D3-D7 background (2.10), 
(2.11), (2.17), and in the Einstein frame with X; = $°, and analyze the stability of 
it by using the linearized equations of motion. Our analysis below mostly follows [33]. 
Here we notice that RR 3-form flux F3 and NSNS 3-form flux H3 vanish on this background, 
and their fluctuations do not mix with that of the other fields in the linear order. Thus we 
start with the action (6.1). The fluctuation of the metric, dilaton, axion, and RR 5-form flux 
dC, = Fs along with their decompositions in terms of spherical harmonics on 9% are summarized 
in (6.2), (6.4), and (6.5). Here we consider the fluctuations satisfying the gauge fixing conditions 
Vehap) = V°hyq = 0 for the metric and V“CargsgL = 0 for the RR 5-form flux F5 = dC4. For 
simplicity, we denote the background metric, dilaton, axion, and 5-form flux by gyn, ¢, x, and 


Fs, respectively. In this section we define @ = a/R° = 4/R. 


A.1 Some Conventions 


Before writing down linearized equations of motion, we summarize some of our conventions. 
We normalized the e tensor on the five-dimensional Lorentzian part Ms and the S° part in 


(2.18) by 


1 
= 01234 __ 
€01234 = yV —9Ms; €E =e 


=9M; 
1 
es6789 = V9ss, ©? = —. (A.1) 
gs5 
In this convention, 
Eniuousiap, e 7 ESRA = —5! ; ETE o a i = 5!. (A.2) 
We also define the ten-dimensional € tensor by 
€0123456789 = V —910 = €01234 ` €56789- (A.3) 
It is also useful to define the Laplacian for Ms and S° by 
g= GV Vos y = 9°? Vag. (A.4) 


A.2 Spherical Harmonics on S° 


In this appendix, we decompose the linearized equations of motion for the fluctuations by using 
the spherical harmonics on $°. Thus we define Y/, ył j Ya B) and Ya ay which represent scalar, 
vector,traceless symmetric, and antisymmetric spherical harmonics, respectively. They satisfy 
the transverse conditions 


ayl _ way —vayl _ 
Vi av ev a ee (A.5) 
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For the vector spherical harmonics, we define the Hodge-de Rham operator A, by AYZ 6 = 


-— RIY. We can define the Hodge-de Rham operators for the other spherical harmonics 


in e similar manner and the eigenvalues of A, on the S% with the radius R are given b 19 ga [32] 


AY" = DY" = ey (k =0,1,2,+*-) (A.6) 
Ayes 1-7) aA (k =1,2,---) (A.7) 
AyYoag) ( y =) (ab) = = i (k = 2,3,4, +) (A.8) 
AyYiag) = ( Y 5) Yaa = U (k=1,2,--). (A.9) 


A.3 Five-form Flux Equation 


Let us first consider the perturbation for the RR. 5-form flux Fs. It satisfies the self-duality 
constraint 


1 
FpaQrst = aroni HE ABCDE. (A.10) 


from which the equation of motion d* F5 = 0 follows automatically. By denoting the fluctuation 
of the 5-form around the background as ôF5 = fs, (A.10) in the linear order is written as 


1 
ÎPQRST = T E fFABCDE +> 3.5] zahm I eporst O” FaBone 
1 
BCD AA 
~ FI PORSTAL E Fa gepeh™®. (A.11) 


By substituting (6.2) and into this equation and then decomposing it by the spherical 


harmonics on S$ 5 we obtain five nn 


or +b! »)| y! =0, (A.12) 


I QI 
SV 11 Os yzpaus T Empzusnans ($r = 


[4V m blousana F Em uonan (10y eB, )| Yd = 0, (A.13) 
iuu + Er pos Vvb al Vee SU, (A.14) 
Viata ia [a] 5 eua One =e est = 0, (A.15) 
Brapa Fna Vorti] Vip’ 4 = (A.16) 


Here & = a/R’ = 4/R. Now, we can simply solve (A.14) and (A.16) algebraically, assuming 
k > 1. Then, we obtain the following three equations from (A.12) and (A.14), (A.13) and (A.16), 


16 Please distinguish the total angular momentum k of the spherical harmonics of S° from the number k of 
D7-branes. 
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and (A.15) respectively: 


[ sta Zh’ — on =| Y! =0, (A.17) 

[Osb] — V’V by + Aybi, — &BI] YJ = 0, (A.18) 
I a YyzV2 pI / J 

[37s Bhans F 2 Em pops biis -A | Yag = 0. (A.19) 


A.4 Einstein Equations 


Let us next consider the perturbation for the Einstein equation. The Einstein equation itself 
can be obtained from the action (6.1) as 


1 1 1 
Run — IMN (z = ze exo x = zoPoo” b — T TaFronsrF?#s" ) 
1 1 1 
-3% OuxOnx = 30m oOne — L I g EMPoRsFy PQRS _ (A.20) 
By using 
1 1 
R- 5e OPXO x — zoo" $ =0. (A.21) 
and 
FporsrFYOPST = 0, (A.22) 


derived from the trace part of (A.20) and the self-duality condition (A.10) respectively, we can 
reduce the Einstein equation (A.20) to a simpler form: 


if 1 
Run — 50 OMXONX — ZðnóðNo a a F MPQRSFy PERS L 0, (A.23) 


By linearly perturbing (A.23) and decomposing it in terms of the spherical harmonics on S5, we 
obtain some equations for the fluctuations around the background. We summarize the resulting 
equations below. 


A.4.1 aß Components 


From (a, 3)-components of (A.23), we obtain the following four equations: 


2 a 
| z +O — = Ya) = 9, (A.24) 
VOB VGYa =0, (A.25) 
v + Z ViaVay¥! = 0, (A.26) 
1 @ 4a i, a orl er 
5 rT + 95" y + io” y + z’ yo eo y = 0. (A.27) 
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A.4.2 pa Components 


In a similar way, from (u, a)-component of (A.23), we obtain the following equations: 


1 V 
z 2Bl — V”V„B}) + IBA, + Sold, + ST ars YJ =0. (A.28) 


1 I 1 I 2 I a oo I 
-5v + 5 ult + ra + zyr? + 96 Se PRP 


1 1 
+ô wz Ow ni + burs o g! |VaY! =0. (A.29) 


A.4.3 pv Components 


Finally, from (u, 7)-component of (A.23), we obtain the following equations: 


1 1 
0 = —-(O0r +0) iw — zVuVo(h +r) + zVuV’ hov + ViVP Ron) + Rupsvh?? 
e? 
+5 (Rahvp ga igh = e°? (OunxXOvX)Y =, -y (AuxOn + OvXOn1) 
1 A x2 x2 
E Pvp — Loboo- -r „eP1P2P3P4P5 7 bI yl = Ipvh + lw 


P1 “ P2P3P4P5 
(A.30) 


A.5 Dilaton and Axion Equations 


Let us then move to the linear perturbation for the dilaton and axion equations of motion. They 


are given by 


Ou(V—g9 On) — V Zget gN OuxOnx = 0, (A.31) 
Ou(V—ger? gg! Anx) = 0. (A.32) 


respectively. When fluctuating around our background, these equations are rewritten as 


(Cn + Dp) + 398-98; (h + 1) — (OrG)(Vuhl") — (Y187)h” 

—2pe?? g?™ (DuxX) + PPR (Owx)? — 267g” (Owx) (Own) =0, (A33) 
(Cr + Oy) + 28-6) (3I + 29" (ue) (wx) + Zdu lh + 2) Bux 
(V h) (uX) — RI (10x) — 2h" (8-4)(Bwx) =0. (A.34) 


A.6 Spectrum for Decoupled Modes 


As we derived the linearized equations of motion for the fluctuations in the appendix [A3] 
and [A.5| we then analyze the spectrum for them. We start with the one for those modes which 
do not mix with other modes in a complicated way. The analysis turns out to be essentially the 
same as the case of AdS; x S° [32]. 
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A.6.1 ¢! Mode from hiap) 


First we consider the scalar mode hap) = gi Yo B) from the fluctuation of the metric. It obeys 
the equation of motion (A.24) and, by using (A.8), we obtain 


( -- OP) o (k = 2,3,4,- --) (A.35) 


R2 


Thus we find its mass m? = EEA, which obviously satisfies the stability condition. 


A.6.2 bi Mode from Crag 


For bl, from the fluctuation of C vag, (A.19) leads to 


(Maxs + Ay)bj,, Yaa) = 0. (A.36) 


The Maxwell operator Max, on Ms is defined by Max,b!, = O,b,, — V”V pb, for the vector 
bl, and we can generalize the definition for teonsors. Since the mass for this mode is given by 


m? =—A, = a (k = 1,2,---), this mode turns out to be stable. 


A.6.3 bl, and Bi, Modes from gua and Cyog-+ 


Let us next consider bi and B3, from the fluctuation of gaa and Cpapy respectively. From (A.18), 
(A.28), and (A.16), we obtain the equations for these modes as 


((Max, + Ay)b/, — aBl)YJ =0, (A.37) 
ã 
((Max, + Ay) Bi — 5 (Max — A), = 0. (A.38) 


or, by denoting in a different expression, as 


bI A =, I 
Mase: ( DEEN "aJ (=o (A.39) 
B, Ay Ay- > B, 


The eigenvalues of 2 x 2 matrix in the second term are 


—~mi=A,-—+1/— - A, . (A.40) 


More explicitly, by using (A.7), we get the masses for these vector fields 


k? -1 (k + 3)(k +5) 
m = -5 3 m —— > (k = 25 3, oe -). (A.41) 


Therefore we can find that these modes are stable, too. 
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A.7 Spectrum for Mixed Modes: h/, nt, and b 


Let us determine the spectrum for the scalar perturbation h!, 7’, and b? here. We first assume 


k > 2 and then find h? = =q! from (A.26). Thus we can rewrite (A.17) and (A.27) as follows: 


b! i —3a b! 
A me E = 


The matrix in the second term is diagonalized and the eigenvalues are given by 


—- m? = Oy — & + ,/a4 — 4420,, (A.43) 
or, more explicitly, by 
k? — Ak 
k+4)(k+8 
„o BEDR (nas 


Even for k = 0,1, we can see that the expression (A.45) is correct. Thus we find the lowest mass 


in this mode is m? = fe when k = 2. This violates the stability bound. 


A.8 Spectrum for Other Modes: h, p, and 7 


Now, the remaining modes are h,,, Y, and 7. Since the analysis of massive modes looks highly 
complicated, below we would like to consider only zero modes on S° of hw and y,n. They 
are useful to the calculations of the viscosity in section [4] For this reason, we generalize the 
background to the black brane metric and write down the linearized equations for the 
fluctuation around it. By taking u — 0, we can reproduce those for the background metric 
(2.18). 


A.8.1 Zero Modes on S° 


We concentrate on the zero modes on S5 i.e. k = 0 modes of the spherical harmonics. Notice 
that in this case, we have b = 0 and by, poy3y4 1S expressed in terms of h and 7 as 
2 
5V yı buonsuans = Re vaHs Mais (h = T) ’ (A.46) 


from (A.10). From (A.27), the mode m satisfies 


( = =) pau (A.47) 


It is also useful to define the Weyl shifted metric H,,, by 


1 
Hw = hw I gIwT, (A.48) 
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which corresponds to the metric perturbation around the five-dimensional background obtained 


by the reduction of S°. Then the Einstein equation (A.30) can be rewritten as follows: 


: (O,d0,~ + O,d0,p) = 0. (A.49) 


4 e2? 
SRG) + =H — e (3 xax) — -y (nxun + ðv xun) — 5 


R2 


Here SRG) is the perturbation of the purely five dimensional Ricci tensor (neglecting the S” 


contributions) due to the metric perturbation H,,,. Notice that when u,v 4 w,r, we obtain the 
simple Einstein equation SRG) + gH yw = 0. This simplification is applicable, for example, to 
the shear modes Hy, and Hy. 


For the dilaton equation of motion (A.33) with k = 0, by using the results 


2 2 H 5u \ 
Tas = -i (1 iis) (1+ sae) = 
2 2 H \— 
Vzrôro = VyOyh = 3" (1 = a} = fe 
4 4 e \ 
Vwowd = g (1 = Hi) = fw, 


_ Hp je NP 
Vid = srg (1- sig) E fu (A.50) 


and the Weyl shifted metric H,,,, we can rewrite it as follows: 


1 
oe + 5-9" (3-H) — fiH” — fia + HY) =F ee oe H" 


1 2 
tar fe +29% fa 4g fy +g” fr) = ave 


1 
+8 eI = Ogos gto) — 2667? gY Ann = 0, (A.51) 


where H = Hii =h+ an. 
On the other hand, for the axion equation of motion (A.34) with k = 0, by using the values 


2k 


vonr- Von 
X X= 3Er 


(A.52) 


it can be rewritten as follows 


4 H 
a + 39 Orn + Bg?” Ow (z + =) — BV ph = 0. (A.53) 


We can take the radial gauge H,,, = 0 and have 5 physical modes for H,,,. Then the dilaton 
and axion equations of motion (A.51)) and (4.53) become 


1 
op + <9" (0-H) — feH" — fa(H™ + H”) 


1 
—fwH?” + gto" fe + 29° fet ge’ fo +9" fr) 


1 
+702 (HY — 2g — zrg) — 28e? getun =0,  (A.54) 
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and 


4, H 
on + zg Orn + Bg” Ow (2¢ + >) — BV ,HY = 0. (A.55) 


In summary, we need to solve the zero mode equations of motion (A.46), (A.47), (4.49), 
(A.54), and (A.55) to find the variables buvpo, T, Y, n, and Huv. 


B Linearized Equations for the Shear and Sound Channel 


In this appendix we summarize the linearized equations for the shear and sound channel. They 


are useful to derive the differential equations for the ae invariant combinations in — 


For this purpose, we consider the equations (A.46| , and (A.55)) derived 


in the appendix [A.8.1] Then we substitute the BINN e to the shear channel 
(4.1) or sound channel (4.3) into them and then derive explicit expressions for the linearized 


equations of motion. We summarize the resulting equations below. 


B.1 Shear Channel 


By considering the fluctuations Hy, and H,, which correspond to the shear channel in the radial 


gauge H,,, = 0, we obtain the equations for them 


H sie )\o. = 32 1 (gfu wHsy) = 0, (B.1) 
a 
Bot ee My Sat + wHzy) = 0, (B.2) 
c2 
qH',, +w aH = 0, (B.3) 


from (t, x)-, (x, y)-, and (x,r)-component of (A.49}), respectively. We can directly show that two 


of the three equations are independent. 


B.2 Sound Channel 


Let us next consider the sound mode. The corresponding fluctuations are Hy, H, a Aye, yy, How, 9 
and fmuzpsuaug ÌN the radial gauge. In this case the linearized equations (A.47), (A.49), and 
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(A.54) are summarized as 


Hy +n! (ef) Hh — In! (ce) (Hy + Fi) 


1 fw? 2 2 
-Z ( z (Hyy + Hii) + au + tte) = 0, 
n 1 (28/3 1 1 1 1 1 1 fa? q? 
Hir gy ln Ge ) Hoz -ln (Cx) (Hi ~~ Hy ~~ Hi,) F 2 E = 2 Hex 
t t x 


Hog Ge Ae Ge Ay A= Ae 


1 
nl (ches?) Hy in el! (= SH) = 0, 


HY + In'(cl#3) Ht, + Sy =0, 


2 2 
CECE 


(Hi, + Hi) + 1H, + In! (cg¢;!)(Hyy + Hi) 


(B.4) 


(B.5) 


1 
4- In’ (c2) Hey + In(c2/) (v — T =0, (B.10) 
w 


Ww (or n 1 
(Ha — Hii) + oe — In! (cre) Hi — In’ (c2/*) (v = T =0, 
t 
" 17/28/37 11/3 1 1 1 1 w2 g 
p +n (cy y — ln (cx Ha — Hy — Hu) + a e e p 
t t T 


27r28/3 (r 1/3 — p) TS 


1 
1 fu? ¢ 32 
a ln (2A?) +> (5 = le Hr O, 
( t <x ) 2 e c2 e 


44 H 11(u3 — 2p?r"!/3 + earl) 
9c? z= 


where Hj; = Hre + Hww. 


C D3-D5 Scaling Solution with F-string Sources 


(B.11) 


(B.12) 


Here we briefly review the D3-D5 scaling solution in type IIB supergravity with F-string sources 


35|. This solution can be regarded as a back-reacted solution dual to a homogeneous baryon 


condensation in M = 4 super Yang-Mills. Each baryon (= D5-brane) situated near the horizon 


r = 0 carries N F-strings which extend into the boundary r = œo due to the string creation [34]. 


Below, we work in the Einstein frame and the supergravity action is obtained by rewriting the 


one in the string frame (ZI) by the Einstein frame metric Gk, = e7%/ “ee 5 
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For this brane setup, it is appropriate to use the following ansatz for the metric 


3 
ds? = —e") A(r)dt? + &b (La) +e") (Al (r)dr? + r7ds%.), (C.1) 
i=1 
where Xs represents a Einstein manifold with the same Ricci curvature as the unit radius S°. 


For fluxes, the ansatz is 


F; = ežblr) . h(r) - dx, A dz ^ dz3, 
Fs = r5e30) f(r)(Qx, + *2x,), (C.2) 


where Qy, is the volume form on X5. We also assume that the dilaton is dependent on the 
radial coordinate r only. If we treat this solution within supergravity, we will have a tadpole 
of the NSNS 3-form flux H3 as is clear from the equation of motion in the presence of the 
Chern-Simons term 


d(e~? x H3) = F3 A Fs, (C.3) 


where F3 and Fs are sourced by the D5 and D3-branes, respectively. Thus we cannot construct 
consistent supergravity solutions under the assumption of symmetry of spacetime implied by 
(C.1). 

To resolve this problem, we notice that the effective F-string charges are generated from the 


D3-brane and D5-brane charges via the Chern-Simons term 
1 
Sos = -qy | Ca MHs A Ps (C.4) 
AK {0 


Therefore, to resolve this problem, we treat the created F-strings as the probe action and add it 
to the supergravity action. The probe action for a single F-string in the Einstein frame is given 
by 

oe -> f drdoe? -Gpr + = I drdo B,- X"3 X”. (C.5) 


By identifying 7 = t and o = r and taking the sum over infinitely many F-strings, we obtain 


the probe action for the created F-strings 
> Sstring(i) = J derdezdzs9x, P Sstrings (C.6) 
i 
where 7 labels i-th F-string and the density of F-strings p is assumed to be constant. In the 


following, we solve the equations of motion derived from the supergravity action with this probe 
action added. 
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C.1 Equation of Motions 
The Bianchi identity for F3 and the equation of motion for F5 are written as 
A,(r°e2") f(r)) =0, O,(e2 h(r)) = 0. (C.7) 


Here we notice that the equation of motion for F3 is automatically satisfied for the ansatz 


introduced above. From these equations, we can define the constants F and H as follows: 
fm) =F. p be 3") h(r)=H- e387), (C.8) 


Then, from (C.3), the density of F-strings turns out to be 


nr H 
2k20 ` 


p= (C.9) 


To derive the dilaton equation of motion, we notice that the F-string action now looks like 


1 1 ec .u 
S'F_string — -5 f atare? V -GE = Dae f etaret #4. (C.10) 


Thus the equation of motion becomes 


Ə, (o (r) Alre O+ $+ P25) = ži e PHPH 4 ESAO)? (0.11) 
We can also derive the Einstein equations for the type IIB supergravity with the F-string 


action. Combined with (C.11), we can summarize the equations of motion as follows: 


A(log A + u)'r5ežbt2e+zuy _ “FHeltictie + ~ beter horde 4 Pree, 


(C.12) 

bf Ades ote guy! = -F Helter T PA ~ = ebtio Hotty (C.13) 
(c+ 2log ry Arnet er = —{FHelt ete B Ege 

+EH ried tei + 8rõe?ct3b+zu (C.14) 

[(log A + u — b — 2¢)' Ar®esbt2e+34y = 0, (C.15) 

2(¢')? +30)? + 108 — 30'c! + 6b" + 10c” — 3b'u! — 5c'u' o% =0. (C.16) 


It is also useful to derive the following equation from a linear combination of (C.12), (C.13), 
and (C.14) 
[log A +u +b + 2c + 4log r) Arex t2e+ au)! = 16r3eret got ae, (C.17) 


In the above discussion, we have derived five equation of motion for five variables A(r), 
u(r), b(r), c(r), o(r). However, we can eliminate one of them, say u(r), by the diffeomorphism 


r—+f=*f(r). In this sense, the independent degree of freedom under the symmetry ansatz is 
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the four variables. Thus we should show the five equations of motion are not over constrained. 
Indeed we can show the following identity from the four equations of motion (C.12), (C.13), 
(C.14), and 

[Acre Heut) S 2 0, (C.18) 


where Eeonstraint(r) is the left-hand side of (C.16). This guarantees that the fifth equation of 


motion Evonstraint(’) = 0 is satisfied if it is vanishing at any particular value of r. 
C.2 Scaling Solutions 
Let us assume the following simple scaling profile for the unknown functions: 
u(r) =urlogr+uo, b(r)= bılogr + bo, c(r)=calogr+o, ¢(r)=¢ilogr+¢o. (C.19) 


At first we also assume A(r) = 1, which corresponds to the extremal case. 
Using the three equations of motion (C.14), (C.15), and (C.16), we find that the following 


coefficients u,,61,c; and ¢, satisfy equations of motion 


2 2 
ui = nj? by = a Ci = —2, Qı = TE (C.20) 


By substituting (C.19) with (C.20) into (C.12), (C.13) and (C.14), we can find four solutions 
for the pair (F, H) in terms of uo, bo, co and ġo as follows: 


83 — v 1081 


F = +e. 


Ho = 2s : s 
m e F, (C.21) 
and 
P= tere V+ 
V3 ? 
g = _MIBITS go-go- p (C.22) 
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Thus, in order to realize the fluxes with arbitrary values and sign, we have only to choose one 
of the four solutions and tune cg and ġo appropriately. We set ug = bo = 0 below. 

As we constructed the extremal solution, generalization to a black brane solution is straight- 
forward by considering more general functions A(r). Let us again assume the ansatz with 
the values (C.20). Then it immediately follows that all equation of motion are satisfied only if 


M 


where M is an arbitrary constant related to the ADM mass of the black brane. We also notice 


that the other profiles are the same as the extremal solutions. 
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In summary we obtain the following metric in the Einstein frame and the dilaton 
2 


dr 
2 _ pu 2 | pbigz2 1 co 2 d(r)-d0 — p% 
ds r4 A(r)dt* + r7'dz* +e (= Cy + (dQs5) ) , œe pe (C.24) 


1 
and, after the redefinition of radial coordinate p = rv, we reach the expression (8.1). 
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